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The phase structure of two-flavor QCD is investigated at real isospin and imaginary quark chemical po- 
tentials by using the Polyakov-loop extended Nambu-Jona-Lasinio model. In the region, parity symmetry is 
spontaneously broken by the pion superfluidity phase transition, whereas charge-conjugation symmetry is spon- 
taneously violated by the Roberge-Weiss transition. The chiral (deconfinement) crossover at zero isospin and 
quark chemical potentials is a remnant of the parity (charge-conjugation) violation. The interplay between the 
parity and charge-conjugation violations are analyzed, and it is investigated how the interplay is related to the 
correlation between the chiral and deconfinement crossovers at zero isospin and quark chemical potentials. 
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I. INTRODUCTION 

The phase diagram of Quantum Chromodynamics (QCD) is 
the key to understanding not only natural phenomena such as 
compact stars and the early universe but also laboratory exper- 
iments such as relativistic heavy-ion collisions. Quantitative 
calculations of the phase diagram with the first-principle lat- 
tice QCD (LQCD) have the well-known sign problem when 
the baryon chemical potential /ib is real [1]; here — 3/iq 
for the quark-number chemical potential ^q. 

The grand canonical partition function Z{^q) of two-flavor 
QCD can be obtained by 

ZifXq)^ JvUdetA{fiq)e-^'', (1) 

where J VU denotes the path integral with respect to gauge 
variable U, Sg stands for the pure gauge action and 

A{Hq) = D^-/^ + Too - ^q74 (2) 

with the gamma matrix 7^ in the Euclidean notation and the 
quark mass niQ. For simplicity, up and down quarks are as- 
sumed to have the same mass toq. When /iq is real, it is easy 
to verify 

^(Mq)^ = -D^-y^ + niQ - /^q74 = 75Zi(-Mq)75- (3) 

This leads to 

{detZ\(/^q)}* = det{75Z\(-^q)75} 

= dctZ\(-/iq) =^detZ\(Mq), (4) 

and hence det /\(/iq) is complex. This causes the sign prob- 
lem. Several approaches have been proposed so far to circum- 
vent the difficulty; for example, the reweighting method [2J, 
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the Taylor expansion method f3\ and the analytic continua- 
tion from imaginary /iq to real fj,q ||4-tl2] . However, those are 
still far from perfection particularly at /iq/T > 1, where T is 
temperature. 

As an approach complementary to LQCD, we can con- 
sider effective models such as the Nambu-Jona-Lasinio (NJL) 
model iflsl - Esll and the Polyakov-loop extended Nambu- 
Jona-Lasinio (PNJL) model |29-72]. The NJL model de- 
scribes the chiral symmetry breaking, but not the confinement 
mechanism. The PNJL model is extended so as to treat both 
the mechanisms approximately by considering the Polyakov 
loop in addition to the chiral condensate as ingredients of the 
model. 

In the NJL-type models, the input parameters are usually 
determined from the pion mass and the pion decay constant 
at vacuum (/iq = and T = 0). It is then highly nontrivial 
whether the models properly predict the dynamics of QCD 
at finite /iq. This should be tested from QCD. Such a test is 
possible at imaginary /iq and finite isospin chemical potential, 
since LQCD has no sign problem there. 

For imaginary quark chemical potential, /iq = iOT, one 
can verify that 

AiiOT)^ = -Z?,7, + Too + Z0T74 = j^^iiOT)^^. (5) 
This leads to 

{det AiiOT)}* = det{75Zi(i6>r)75} = det A{ieT), (6) 

and hence dot A(i0T) is real. LQCD has thus no sign prob- 
lem at imaginary /iq. 

At imaginary /iq, the QCD partition function has the 
Roberge-Weiss (RW) periodicity, i.e. the periodicity of 27r/3 
in 0. The RW periodicity is a remnant of the Z3 symme- 
try in the pure gauge limit. This periodicity can be reinter- 
preted as the extended Z3 symmetry |47, 48, 52,.54J; see ( [30l l- 
( I33I 1 for the relation between the RW periodicity and the ex- 
tended Z3 symmetry. At higher T, three Z3 vacua emerge 
one by one |73] as increases. This mechanism guaran- 
tees the extended Z3 symmetry (the RW periodicity) at higher 
T. The mechanism induces the first-order phase transition at 
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O = {2k + l)7r/3, where k is an integer fl^. This phase tran- 
sition is called the RW transition. The partition function of the 
PNJL model has the extended Z3 symmetry and hence yields 
the RW periodicity and the RW transition (4§, H |70M72|l. 
Eventually, the PNJL model qualitatively reproduces LQCD 
dataH. 

LQCD simulations at imaginary /iq show that the chiral and 
deconfinement crossovers take place simultaneously [4-12]. 
Such a strong correlation between the transitions is not seen 
in the original PNJL model iHllUl^l- We then extended the 
PNJL model to solve this problem. In the new model called 
the entanglement PNJL (EPNJL) model, the strength of the 
four-quark vertex depends on the Polykov loop ll63ll . In the 
EPNJL model, the chiral and deconfinement crossovers coin- 
cide with each other, so that the EPNJL model yields good 
agreement with LQCD data [6^ . 

The reliability of the PNJL and EPNJL models can be tested 
at finite isospin chemical potential (/ijso), since LQCD has no 
sign problem there 1.74,] . For later convenience, we use the 
"modified" isospin chemical potential /ij = /iiso/2 instead of 
/^iso = A'm — Md, where /i„ and are chemical potentials for 
up and down quarks, respectively. The QCD partition function 
at finite is obtained by 

Z{fii) = JVU dot A{nu) dct A{nd)e^^'' 

= JVUdctA{^ii)detA{-^ii)e-^^. (7) 

When /i„ and /id are pure imaginary, both det A{fiq) and 
det /\(/iq) are real as shown in (|6]l. When /i„ and fid are real, 
it is possible to prove that 

{detA{ni)detA{~ni)}*= det A{ni) det A{- m) . (8) 

The product of the two determinants is thus real, whereas each 
of the determinants is not. LQCD has hence no sign prob- 
lem for both real and imaginary /^i. Actually, LQCD data are 
available there; see Refs. H for real /ij and Refs. iS H 
for imaginary /ij. 

The PNJL model ex plai ns the LQCD data qualitatively not 
only at real fj,i ll36i I37il62ll and but also at imaginary /ii IstII : 
the eight-quark interaction newly added improves the agree- 
ment between the PNJL result and the LQCD data iH]. The 
EPNJL model reproduces the LQCD data quantitatively with 
no free parameter; note that all the parameters in the EPNJL 
model are fixed at vacuum and at imaginary chemical po- 
tential. As for real /ij, there exists a tricritical point (TCP), 
i.e. a meeting point between the first- and second-order pion- 
superfluidity phase-transition lines. The location of the TCP 
is predicted by the EPNJL model |63] as well as the PNJL 
model |36, 62], the chiral perturbation theory Jt^I, the strong 
coupHng QCD |76] and so on IttI]. 

As for real /iq, LQCD data are available at small /iq/T 
with the Taylor expansion method |3]. The EPNJL model re- 
produces the LQCD data quantitatively with no free parame- 
ter [67]. The EPNJL model predicts |63] that there is a critical 
endpoint (CEP) ifTil FtsIi of the first-order chiral transition in 
the fiq-T plane. 

LQCD has no sign problem also at real isospin and imagi- 



nary quark chemical potentials, because 

{det Z\(/ii + i0T) det A{-fii + iOT)}* 
= det A{ni + iOT) det A{-fii + ieT) (9) 

for /iu = /ii + iOT and /i^ — — /ij + iOT. The reliabil- 
ity of the PNJL and EPNJL models can be checked also in 
this case, although no LQCD data is available at the present 
stage. It should be noted that the LQCD calculation with 
real ui and imaginary /iq is equivalent to the LQCD calcu- 
lation with complex quark number chemical potential under 
the phase quenched approximation [79]. 

There is a debate whether the chiral and deconfinement 
transitions at /iq = /ij = coincide or not; see Ref. fsO*] and 
references therein. If the two transitions do not coincide, ex- 
otic phases such as the constituent-quark phase ll8ll[82ll or the 
quarkyonic phase ||46[ Isil HH HH 115] may appear However, 
the chiral and deconfinement transitions are confirmed to be 
crossover with LQCD ll85l[86ll . so that transition temperatures 
of the two crossovers are not well defined. The coincidence 
problem thus includes conceptual difficulty. 

As a way of circumventing the conceptual difficulty in the 
coincidence problem, we should consider exact phase transi- 
tions relevant to the chiral and deconfinement crossovers. This 
is really possible. It is reported with the EPNJL model ll66|] 
that the chiral and deconfinement transitions coincide when 
the parity (P) restoration and the charge-conjugation (C) vi- 
olation occur simultaneously at 6 = = tt, where 6 is the 
parameter of the (^-vacuum ll25l - l28ll . The chiral and decon- 
finement crossovers at /tq = /xi = are remnants of the P 
restoration and the C violation at 9 ~ ~ n, respectively. 
LQCD has, however, the sign problem at finite 6 and hence it 
is difficult to check the validity of the model prediction with 
LQCD. Such a problem does not appear at /tj > and ^ tt, 
as mentioned above. 

In this paper, we investigate the interplay between the pion- 
superfluidity and RW phase transitions in the region of = tt 
and /ii > m„/2, using the PNJL and EPNJL models; here 
is the pion mass in the vacuum. In the region, the P vio- 
lation due to the pion-superfluidity transition occurs at lower 
T, whereas the C violation due to the RW transition takes 
place at higher T. These transitions are exact phase transi- 
tions and hence we can define the order parameters without 
any ambiguity. We also investigate the relation between the 
chiral transition at /ij = 6* = and the P violation at 6* = 
and fii > by varing /ij and the relation between the 

deconfinement transition at /ij = = and the C violation 
at 6* = TT and /ij = by varing 0. Finally we discuss how 
the correlation between the P and C violations at = tt and 
III > TOjr/S is related to the correlation between the chiral 
and deconfinement crossovers at /tj = = 0. In future, we 
can check the model prediction made in this paper by using 
LQCD, since LQCD has no sign problem in the region and 
hence LQCD simulations are feasible. 

This paper is organized as follows. In section II, we explain 
the PNJL and EPNJL models. In section III, numerical results 
are shown. Section IV is devoted to summary. 
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II. PNJL MODEL 

The two-flavor PNJL Lagrangian in Euclidean space-time 

is 

+ U{$[A],$[A]*,T), (10) 

where D^, = iA„ and A„ = S^^/^gA^^a ^ with the gauge 
field a, the Gell-Mann matrix Aa and the gauge coupling 
g. In the NJL sector Gg denotes the coupling constant of 
the scalar-type four-quark interaction. As is seen later, the 
Polyakov potential Uisa. function of the Polyakov loop ^ and 
its Hermitian conjugate <P* . 

The quark mass matrix toq is given by ttiq = 
diag(mo, mo)- The chemical potential matrix ji is defined by 
ji — diag(/i„, Md) with the u-quark (d-quark) number chem- 
ical potential /i„ (/id)- This is equivalent to introducing the 
baryon and isospin chemical potentials, ji-Q and jiiso, coupled 
respectively to the baryon charge B and to the isospin charge 



^J■ = MqTo + /ilT3 



(11) 



with 



Mil + Mb Mu ^ Md Miso 

Mq = = W = = — ' (12) 

where tq is the unit matrix and (i = 1, 2, 3) are the Pauli 
matrices in the flavor space. Note that /iq is the quark chem- 
ical potential and mi is half the isospin chemical potential 
(Miso)- In the limit of mo — iii = 0, the PNJL Lagrangian 
has the SUl{2) x SUr{2) x U^{1) x S';7c(3) symmetry. For 
mo 7^ Oand/ii ^ 0, it is reduced to C/13 (1) x {7v(l) x S'[/c(3). 
The Polyakov loop operator and its Hermitian conjugate 
are defined as 



1 

iV 



TrL, <pt 



1 

N 



TrLt 



with 



i(x) = V exp 



(iT^4(x, r) 



(13) 



(14) 



where V is the path ordering and = iAq. In the PNJL 
model, the vacuum expectation values, (p ~ (3) and <P* — 
{'P^), are treated as classical variables. In the Polyakov gauge, 
L can be written in a diagonal form in color space ll3Tll : 

^ = g»/?(03A3+08As) ^ diag(e*^'^",e'^'^',e*^'^<=), (15) 

where = 03 + ^s/VS, "/"b = -03 + 4>8/Vi and 0c = 
-(0a+0&) = -208/V3. 

The Polyakov loop ^ is an exact order parameter of the 
spontaneous Z3 symmetry breaking in the pure gauge theory. 
Although the Z3 symmetry is not exact in the system with 
dynamical quarks, it still seems to be a good indicator of the 



deconfinement phase transition. Therefore, we use 'P to define 
the deconfinement phase transition. 

The spontaneous breakings of the chiral and the Ui,j{l) 
symmetry are described by the chiral condensate a = (qq) 
and the charged pion condensate |.2L.22..36„,62.1 



(16) 



where t± — (ti ± ir2)/V2- Since the phase (p represents the 
direction of the C/i., (1) symmetry breaking, we take ip = for 
convenience. The pion condensate is then expressed by 



vr = {qi-f5riq). 

Making the mean field (MF) approximation 
can obtain the MF Lagrangian as 



(17) 



with 



M 
N 



niQ - 2GsCr, 

-2GsTT. 



U 



(18) 



(19) 
(20) 



Performing the path integral in the PNJL partition function 



^PNJL 



DqDq exp 



(TxC 



MF 



(21) 



we can get the thermodynamic potential Q (per unit volume). 



-Tln(ZpNjL)/V^--2^ / 

i=± 

\ In [1 + 3(<? + <?*e-/5^^r(p))e-/^^^r(p) 



3£;,(p) 



-3/3S, (p)l 



1 



In [1 + ii'P* + <Pe-~PEt{v)y-PEt{v) + g-3/3B+(p)] 



U 



with i?±(p) = -B±(p) ± Mq' where 



(22) 



(23) 



for i?(p) = yP^ + M^. On the right-hand side of (l22T i. only 
the first term diverges, and it is then regularized by the three- 
dimensional momentum cutoff A 1,31. .34)] . 

We use U of Ref. iH fliat is fitted to LQCD data in the 
pure gauge flieory at finite T iIstIIssII : 



a(T) 



+ b{T) ln(l - Qm* + 4(<p3 + <|>*3) _ 3(^<p*)2) 



a{T) = flo + ai + a2 



rn\2 /To\3 
6(r)=63(^' 



(24) 

3 

(25) 
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where parameters are summarized in Table J] The Polyakov 
potential yields a first-order deconfinement phase transition at 
T = To in the pure gauge theory. The original value of Tq is 
270 Me V determined from the pure gauge LQCD data, but the 
PNJL model with this value of Tq yields a larger value of the 
pseudocritical temperature Tc at zero chemical potential than 
Tc = 173±8 MeV predicted by the full LQCD simulation fs?, 
IMIM]. We then reset Tq to 212 MeV {SM so as to reproduce 
the LQCD result. 



ao 


ai 02 


b:i 


3.51 


-2.47 15.2 


-1.75 



TABLE I: Summary of the parameter set in the Polyakov-potential 
sector determined in Ref. 13511 . All parameters are dimensionless. 



Table shows parameters in the NJL sector used in the 
present analyses. This set can reproduce the pion decay con- 
stant = 93.3 MeV and the pion mass to^ = 138 MeV at 
vacuum (T = iiq = ~ 0). 



Gs 



mo 



A 



5.498 [GeV~^] 5.5 [MeV] 631.5 [MeV] 



TABLE II: Summary of parameters in the NJL sector. 



The classical variables X = 
by the stationary conditions 

df2/dX = 0. 



a and tt are determined 



(26) 



The solutions to the stationary conditions do not give the 
global minimum of i7 necessarily. There is a possibility that 
they yield a local minimum or even a maximum. We then have 
checked that the solutions yield the global minimum when the 
solutions X{T, /iq, /ii) are inserted into (l22b . 

The thermodynamic potential (|22] | is invariant under the 
parity transformation 



TT — )• — TT. 



(27) 



Parity odd quantity such as the pion condensate is an order pa- 
rameter of P symmetry. When pion condensates, P symmetry 
is spontaneously broken |74]. 

When Hi ^ 0, charge-co njug ation symmetry is explicitly 
broken even at /iq = fssi l66ll . since u and d quarks have 
different electric charge. In the present system where the elec- 
tromagnetic field is switched off, however, we can neglect the 
electric charges and consider only a charge-conjugation trans- 
formation associated with baryon and color charges. We re- 
fer to symmetry under the transformation as partial charge- 
conjugation symmetry (C) in this paper. As far as fl of ( |22] | 
is concerned, the transformation is simply represented by 

0^-0 (28) 

with (\) the phase of the Polyakov loop ^. When /iq = iOT = 
in-KT for any integer n, fl is C-invariant. As shown below, C 



symmetry is spontaneously broken at higher temperature for 
odd n. 

The thermodynamic potential fl of (l22b has a trivial pe- 
riodicity of 27r in O. We then mainly consider one circle, 
< < 27r. In addition, fl has the RW periodicity iItI . 
The thermodynamical potential fl is not invariant under the 
Z3 transformation. 



(29) 



for any integer fc, while U of ( l24b is invariant. Instead of the 
Z3 symmetry, however, J? is invariant under the extended Z3 
transformation ll48ll. 



(30) 



The extended Z3 symmetry guarantees the RW periodicity, 
as shown below. The thermodynamic potential (l22b can be 
rewritten with the modified Polyakov loop 



If' = <Pe 



10 



(31) 



invariant under the extended-Z3 transformation: 



=± ■ 



+ iln [1 + 3i?'e-'3^'(P) 

P 

+ ^ln [1 + 3<f *e-''^'(P) 



rp4 



+ b{T) ln(l - m^* + 4:{<P^e-^'^ + <F*^e^'^) - 3{^^*f) 

(32) 

Obviously, J7 is extended-Zs invariant, since it depends on 
only through the factor e'^*®. The factor also guarantees that 



(33) 



Thus, fl{0) has a shorter period 27r/3 in 0. This periodic- 
ity was first found by Roberge and Weiss [73] for QCD and 
now it is called the Roberge-Weiss (RW) periodicity. Using 
the perturbation theory and the strong-coupling lattice theory, 
Roberge and Weiss also found that a first-order phase transi- 
tion occurs at 6* = 7r/3 mod 27r/3 when T is higher than some 
critical value Tq. The transition is also called the Roberge- 
Weiss transition. After the theoretical prediction, the RW pe- 
riodicity and the RW transition were confirmed by LQCD JJ- 
12]. Recently, the RW periodicity and the RW transition were 
also confirmed by Holographic QCD |91]. The RW transition 
induces the C breaking Issl 1^1 for /ij = and the C break- 
ing for Hi 7^ 0. The C and C breakings occur at = tt 
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and the Z3 images of the breakings also appear at = m/i 
and z57r/3. It is convenient to use the phase ijj of the modi- 
fied Polyakov loop ^ as an order parameter of the C and C 
breakings |[55il66ll . When C or C symmetry is preserved, V' 
vanishes; here we consider ij} in the region — tt < ^ < tt. 

In LQCD simulations at imaginary /iq, the chiral and de- 
confinement transitions take place simultaneously, although 
they are crossover [4- 121. Such a strong correlation can not 
be reproduced by the PNJL model[48i jsl El- This prob- 
lem is solved by the entanglement PNJL (EPNJL) model ll63ll . 
In the EPNJL model, the four-quark vertex depends on the 
Polykov loop. In QCD, the four-quark vertex is generated by 
the gluon propagation between two quarks. If the gluon has 
a finite vacuum expectation value in its time component, the 
gluon propagation can depend on ^ through the vacuum ex- 
pectation value [63]. In fact, recent calculations ll^| -|94ll of 
the exact renormalization group equation (ERGE) ll95ll sug- 
gest that mixing interactions between a and are induced. It 
is highly expected that the functional form and the strength of 
the entanglement vertex are determined in future by ERGE. 
Ref. f65], the following (p dependence of G is assumed by re- 
specting the chiral symmetry, P symmetry, C symmetry and o 
the extended Z3 symmetry: ^ 



G{$) = G[l - ai^<P* - a2(<? + <P )] 



(34) 



In the EPNJL model, the chiral and deconfinement crossovers 
coincide at zero chemical potential. The EPNJL model with 
the parameter set, ai — a2 — 0.2 and Tq — 190 MeV, can 
reproduce LQCD data at imaginary chemical potential ll63ll . 
real isospin chemical potential |63] and small /iq/T ll67ll . The 
agreement between LQCD data and the EPNJL result persists 
also under strong magnetic field ll65i l96ll . We then use the 
EPNJL model with this parameter set also in this paper 

The entanglement vertex G(<?) agrees with the parameter 
G at T = since ^ — there. At vacuum the EPNJL model^0H 
is then reduced to the NJL model; note that the PNJL model ^ 
agrees with the NJL model at vacuum. As a merit of this prop- 
erty, the NJL sector of the EPNJL model has the same values 
of parameters as the NJL model. 

Very recently, it was shown that the coincidence problem 
between the chiral and deconfinement transitions does not oc- 
cur also in the nonlocal PNJL model |71]. In the model, the 
entanglement between a and <P naturally appears in the renor- 
malization factor. It is then interesting as a future work to 
study the relation between the nonlocal PNJL and the EPNJL 
model. 



m. NUMERICAL RESULTS 



chiral restoration and the deconfinement transition are seen to 
be crossover. The chiral restoration starts at the same temper- 
ature (T « 170 MeV) as the deconfinement transition does, 
but the former transition is slower than the latter transition. 
Consequently, the pseudocritical temperature of the chi- 
ral restoration becomes higher than that of the deconfi- 
ment transition. Actually, it is found that Tcr = 216 MeV 
and T<i, — 173 MeV |52], when the pseudocritical tempera- 
tures are defined by the peak positions of a and respec- 
tively [45]. Figur^TJb) shows T dependence of a and (p cal- 
culated with the EPNJL model. The chiral restoration occurs 
with the same speed as the deconfinement transition, so that 
= Tg, = 173 MeV [63]. 



o 



0.5 - 








0.1 0.15 0.2 0.25 0.3 

T(GeV) 




Fig. 1: (color online). T dependence of the chiral condensate a and 
the Polyakov loop $ at /ii = and fiq — calculated with (a) 
the PNJL model and (b) the EPNJL model. The solid (dashed) line 
represents a {<!>). Here, a is normalized by the value ao at vacuum. 



A. Results at = and fiq = 

In this subsection, we consider the case of /xi = /Xq = and 
Too = 5.5 MeV. The chiral condensate a and the Polykov loop 
<P are approximate order parameters of the chiral and decon- 
finement transitions, respectively. Figur^H^) shows T depen- 
dence of (7 and <P calculated with the PNJL model. Both the 



When two transitions are first order, they might coincide ex- 
actly iIqtIi . At zero /iq, however, the chiral and deconfinement 
transitions are found to be crossover lf85[f86ll . There is hence 
no a priori reason why the two crossovers coincide exactly. In 
LQCD simulations at zero chemical potential [85, 98, 99], ac- 
tually, there is a debate whether the transitions really coincide 
or not; see Ref. 18Q[ and references therein. 
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The coincidence problem has conceptual difficulty, since 

CT and are approximate order parameters of the chiral and 

deconfinement transitions. This point can be circumvented 

by considering exact phase transitions relevant to the chiral^ 

and deconfinement crossovers. In the previous work ll66ll . we 

found that the chiral and deconfinement transitions coincide— 

o 

when the P restoration and the C violation occur simultane- ^ 
ously at 6* = = vr. The chiral restoration at 6 = = ^ 05- 
is a remnant of the P restoration at 6 ~ O ~ n, whereas the 
deconfinement crossover at = 6* = is a remnant of the C—^ 
violation at 9 = = tt. Note that the P restoration and C ^ 
violation are exact phase transitions. This theoretical predic-"^ 
tion, however, cannot be confirmed by LQCD, since LQCD 
has the sign problem at finite 9. The region of /^i > m7r/2 
and = TT that we consider in this paper does not suffer from 
the sign problem. 



(a) 

\ 

: A 


1 1 




1 1 1 



0.1 



0.15 



0.2 

T(GeV) 



0.25 
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B. Results at finite fii and /i, = 

When > the pion condensation becomes finite 

at low temperature, so that P symmetry is spontaneously bro- 
ken there ll74ll . In contrast, the pion condensate tt is zero af 
Mi ^ "i7r/2, so that P symmetry is preserved there. Figure- 
|2] shows T dependence of the pion condensate at /ii = TOtt 
and fiq = 0. The pion condensation disappears at T > Tp = 
201(170) MeV in the PNJL (EPNJL) model. 

Figure[3]shows the approximate order parameter Vct^ + t? 
of the chiral transition as a function of T and fii calculated 
with the EPNJL model, where — 0. As for /ij > 711^^/2 
where P symmetry is broken at small T and preserved at 
large T, y/a^ + n'^ has either a c usp or a discontinuity, as ex- 
pected. These singularities in Vc^ + tt-^ come from the fact 
that the pion condensate tt has either a cusp or a discontinu- 
ity there; more precisely, the P restoration is second order at 
m^/2 < /ii < 95 MeV and first order at /ii > 95 MeV. 
As for fii < TO7r/2 where the pion condensation vanishes 
and hence P symmetry is preserved, in contrast, y/a^ + n^ 
changes smoothly. However, the change is still rapid, since 
a has a rapid change. The rapid change of a is nothing but 
the chiral restoration at zero and small iij. Thus, the chiral 
restoration at zero fii can be regarded as a remnant of the P 
restoration at fii > m^/2. 



C. Results at fii = and O — n 

When — TT, the RW transition take places at higher T. 
Hence C symmetry is spontaneously broken there |55, 66]. As 
an order parameter of the C violation, we can consider ^-odd 
quantities such as the quark number density Uq or the phase ip 
of the modified Polyakov loop. Figure|4]shows T dependence 
of a, \tp\ and |?iq| at /ii = and = tt. Note that rtq is pure 
imaginary and hence |nq| = |Im(nq)|. In the PNJL (EPNJL) 
model, the C violation occurs at T > Tc = 189(185) MeV, 
since riq and tp are finite there. 

Figure |5] shows T and dependences of the approximate 
order parameter |^| of the deconfinement transition calculated 



o 



o 




Fig. 2: (color online). T dependence of chiral condensate a, pion 
condensation tt and Polyakov loop $ at = and /Xq = calcu- 
lated with (a) the PNJL model and (b) the EPNJL model. The solid 
(dashed, dotted) line represents tt (#, a). Here a and tt are normal- 
ized by IttoI, where ttq = tt{T — 0). 



with the EPNJL model; here the case of /ii = is considered. 
Since is ^-even and has the RW periodicity, we can con- 
sider the region < 6* < 7r/3 only. There appears a first order 
phase transition at 6* = tt/3. This is the RW phase transition 
in which C symmetry is spontaneously broken. The transi- 
tion becomes crossover as decreases from tt/3. However, 
a smooth but rapid change remains even at = 0. The rapid 
change is nothing but the deconfinement transition at zero 0. 
Thus, the deconfinement transition at zero Hq is a rernnant of 
the first-order RW phase transition at 69 = 7r/3 1551, l66ll . 



D. Results at finite ni and 

Figure|6]shows T dependence of the pion condensate tt and 
the phase of If" at /ii = m.^ and — tt. In the PNJL 
(EPNL) model, P symmetry is spontaneously broken below 
Tp — 250(196)MeV as shown by a finite value of the pion 
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Fig. 3: (color online). The order parameter \/ a'^ -\- -k^ of the chiral 
transition as a function of T and plane calculated with the EPNJL 
model. Here we consider the case of /iq = 0. 



condensate. In contrast, C symmetry is spontaneously vi- ^ 
olated above = 188(184)MeV as indicated by a finite 
value of ij}. Thus P and/or C symmetries are always bro- o 
ken at any temperature. This situation is the same as that at.C 
61 = = TT 166]. The difference AT = Tp - is smaller b 
in the EPNJL model than in the PNJL model. The entangle- 
ment interaction thus makes AT very small, although the two 
transitions do not coincide exactly. Note that at /ij ~ TOtt the 
P restoration is second order whereas the C violation is first 
order in both the PNJL and EPNJL models. Comparing Fig.|6] 
with Fig. 121 one can see that Tp goes up as O increases. 

Figure Q shows the same quantities as Fig. |6l except the 
case of III = 600MeV is considered. In both the PNJL and 
EPNL models, the P restoration and the C violation occur at 
the same temperature, more precisely Tp = T^ = 160 MeV 
for the PNJL model and 147 MeV for the EPNJL model. This 
is because the P restoration becomes sharper and eventually 
first-order as increases. Thus the first-order nature of the P 
restoration makes the C violation strong first-order and con- 
sequently attracts the C violation. 

Figure |8] shows the phase diagram in the /ij-T plane at 
= 71. The P violation due to the pion condensation oc- 
cur at /ii > TO7r/2 and lower T. The P restoration with re- 
spect to increasing T is second-order at /ij < 480(450) MeV 
and first-order at /ij > 480(450) MeV in the PNJL (EP- 
NJL) model. The C violation is always first order. When 
/ii > 600(545) MeV, the P restoration coincides with the C 
violation in the PNJL(EPNJL) model, as already mentioned 
in Fig.|2] The region where P and C symmetries are violated 
simultaneously is tiny in the EPNJL model compared with the 
case of the PNJL model. 

Figure |9] shows O dependence of the pion condensation tt 
and the imaginary part of the quark number density Hq cal- 
culated with the PNJL model; note that riq is pure imaginary 
in this case. In panel (a) where T = 180 [MeV] < Tc = 
188 [MeV], both tt and Uq are smooth functions of and 
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Fig. 4: (color online). T dependence of a, \^f}\ and \nq\ at /ii = 
and = TT calculated with (a) the PNJL model and (b) the EPNJL 
model. The solid (dashed, dotted) line represents a (\4>\, l'^?!)- Here, 
a is normalized by the value ao at vacuum and Uq is by NfT^, where 
TV; =2. 



no RW phase transition is seen at 6* = tt/3 mod 27r/3. In 
panel (b) where T = 220 [MeV] > Tc, the 6>-even quan- 
tity TT has a cusp at 6* = 7r/3 mod 2tt/3, whereas the ^-odd 
quantity Uq is discontinuous there |54]. This singularity is 
called the RW phase transition. The pion condensate is fi- 
nite only around = tt/3 mod 27r/3 and vanishes when 
Tp < 220 MeV. This behavior of tt resembles that of a in 
the case of toq = Mi = 114^ where a and tt are completely 
symmetric. 

Figure [TO] shows the phase diagram in the 0-T plane at 
III = rriTr calculated with (a) the PNJL model and (b) the 
EPNJL model. The phase diagram is ©-even and has the RW 
periodicity. The RW transition at 6* = tt/B mod 27r/3 is first 
order. The deconfinement transition is first-order in the vicin- 
ity of the endpoint of the RW transition, crossover for other 
0. The P restoration is always second order in the PNJL 
model. In the EPNJL model, the P restoration is first-order at 
0/{tt/S) = -0.3 - 0.3, 1.7 - 2.3, 3.7 - 4.3 and second 
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Fig. 5: (color online). The approximate order parameter \$\ of the 
deconfinement transition as a function of T and 6 calculated with the 
EPNJL model. Here the case of = is taken. 
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Fig. 6: (color online). T dependence of the pion condensate n and the 
phase i/; of tf' at /ii — m,r and = n calculated with (a) the PNJL 
model and (b) the EPNJL model. The solid (dashed) line represents 
|7r| (\'ip\) and tt is normalized by the value ttq at T = 0. 
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Fig. 7: (color online). T dependence of the pion condensate n and 
the phase -0 of i?' at /ii = 600 MeV and — tt calculated with 
(a) the PNJL model and (b) the EPNJL model. See Fig. |6] for the 
definition of lines. 



order for other 0. Thus the imaginary quark chemical po- 
tential makes the P restoration weaker. In the regions where 
the P restoration is first-order, the P restoration and the de- 
confinement transition almost coincide with each other in the 
EPNJL model. 



IV. SUMMARY 

We have studied the interplay between the pion- 
superfluidity and RW phase transitions at real isospin and 
imaginary quark chemical potential, using the PNJL and EP- 
NJL models. In the region of 6* = tt and /ii > par- 
ity symmetry (P) and charge-conjugation symmetry {C) are 
spontaneously broken by the pion-superfluidity and RW tran- 
sitions, respectively. The chiral and deconfinement crossovers 
at zero isospin and quark chemical potentials are remnants 
of these exact phase transitions, i.e. the P restoration and 
the C violation, respectively. The chiral and deconfine- 
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Fig. 8: (color online). Phase diagram in the fii-T plane at = tt cal- 
culated with with (a) the PNJL model and (b) the EPNJL model. The 
dashed (solid) line represents the second (first) order P transition, 
while the dotted line represents the first order C transition. 



ment crossovers almost coincide when AT ^ Tp — is 
small. The interplay between the chiral and deconfinement 
crossovers at zero isospin and quark chemical potential is thus 
determined by the correlation between the P restoration and 
the C violation at > m„/2 and = tt. 

At /ii > AniTr, the P restoration becomes first-order and 
hence the P restoration and the C violation take place simul- 
taneously. The P restoration and the C violation are strongly 
correlated with each other in this region. The correlation is 
weakened as /ij decreases from 4to^. At /ij ss m-^ ^ Stott, 
the critical temperature of the P restoration is larger than that 
of the C violation by 10 MeV in the case of the EPNJL model 
and by 60 MeV in the case of the PNJL model. As a con- 
sequence of this property, at zero isospin and quark chemical 
potentials, the pseudocritical temperature of the chiral transi- 
tion is larger than that of the deconfinement transition by only 
a few MeV in the case of the EPNJL model and by 50 MeV 
in the case of the PNJL model. The correlation between the 
P restoration and the C violation at /ij « ^ 2m-^ and 



2 4 
0/(71/3) 



Fig. 9: (color online). dependence of the pion condensate tt and 
the imaginary part Im(nq) of quark number density calculated with 
the PNJL model. Here the case of /ii — is taken. Quantities are 
shown in GeV'^. Panel (a) is for the case of T = 180 MeV and panel 
(b) is for the case of T = 220 MeV. The solid (dotted) line represents 
TT (Im(nq)). 



O — TT is thus the key to determining the interplay between 
the chiral restoration and the deconfinement transition at zero 
isospin and quark chemical potentials. 

The theoretical prediction at 6> = tt and finite /ii can be 
checked by LQCD, since LQCD has no sign problem there. 
Furthermore, the P restoration and the C violation are exact 
phase transitions and hence we can define the order param- 
eters exactly. In future, we can thus determine the critical 
temperatures of the exact phase transitions without any am- 
biguity by using LQCD. LQCD analyses along this line may 
give an insight into the coincidence problem between the chi- 
ral restoration and the deconfinement transition at zero isospin 
and quark number chemical potentials. Furthermore, we may 
get knowledge of the phase diagram and the equation of state 
at real /ij and real /iq by the extrapolation from real /ii and 
imaginary /iq. 

It may be also interesting to study the QCD phase diagram 
in the region where all of 0, /ij, 9 and T are finite. For 
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Fig. 10: (color online). Phase diagram in the 9-T plane at /ii = m-n- 
calculated with (a) the PNJL model and (b) the EPNJL model. The 
dot-dashed lines stand for the RW transition, and the dotted lines rep- 
resent the deconfinement transition. The solid (dashed) lines denote 
the P restoration of first-order (second-order). 



= /ii = = and small or zero T, only a is finite. For 
large jii and low T, vr is the largest condensate instead of a. At 
large 6 and low T, meanwhile, rj is the largest one. For large 
values of /ij and 9, it was shown by the NJL model |26] that, 
instead of tt or 77, the charged-scalar and isovector condensate 
can condense at low T. In this case, the Ui^{l) symme- 
try under the transformation induced by the isospin generator 
T3 is spontaneously broken, whereas parity symmetry is not 
broken. In principle, one can discuss the coincidence between 
the Ui^{l) breaking and the RW transition by using the PNJL 
model. The calculation, however, is quite time-consuming. 
This analysis is then an interesting future work. 
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